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f~^ ' I review the recent progress in solution of the evolution equations of the three particle hadron 

0\J , distribution amplitudes. 



The three particle operators appear in QCD in the studies of the twist - 3 parton distribu- 



QQ , tionsBB'u and tedal ~ 3 meson wave functions QU, and more interestingly, leading twist nucleon 

'nJ" ' wave functions Bffl. For example, the Q^ dependence of a baryon wave function is given to one 



loop accuracy by the following expression 



rrs ' / a (O) \'^^'i/^'^ 

■^ ■ (PiXi, Q^) = a:iX2X3 ^ (t)N,q PN,q{xi) ( ' . ) , (1) 

^' N,q \as[Qo)J 

o ■ where bo = 11/3 Nc — 2/3nf. The anomalous dimensions of the three quark operators, ^N,q, and 

D . polynomials, PN,q{xi) arise from the solution of the eigenvalue problem for a certain integral 

operator (Hamiltonian) H, which effectively describes the dynamics of the three particle system 

with a pairwise interaction. It has been observed recently H that for the particular twist-3 

^ . operators the corresponding three particle systems are intrisically related to the Heisenberg spin 

H \ magnet, which is known to be integrable quantum mechanical model and their solution can be 

found by applying powerful method of Integrable Models. 

It is well known that the conformal symrnetry of the QCD Lagrangian reveals itself as the 
SL{2, R) invariance of the evolutions equationslly'td. As a consequence the pairwise Hamiltonians 
depend on the corresponding two - particle Casimir operators only. 



H — 2_^ Hik, Hik — h{Jik), Jik{Jik + 1) — Lik — i^i + ^k) , (2) 

ik 

where Li is the SL{2, R) generators associated with i-th particle. Then the explicit calculationscJ 
yield that up to unessential constant the Hamiltonian governing the Q^ evolution of the wave 
function of baryon with the maximal helicity A = 3/2 

H3/2 = 2 (1 + 1^) T.i^{J^k + 1) - V'(2)] + ^Cf (3) 

^ "^^ i<k 

coincides with the Hamiltonian of the noncompact XXXs=-i Heisenberg spin magnet. This 
Hamiltonian possesses an additional integral of motion (conserved charge) 

Q = 2 [-^12' ^23] = 'i'^did2d3Zi2Z23Z31 , (4) 
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Figure 1: The spectrum of eigenvalues for the Hamiltonian Tti2. The lines of the largest and the smallest 
eigenvalues of 7^32 are indicated by dots for comparison. 



that commutes with -^3/2 and with SL{2) generators. It is interesting to note that this model 
has already been encountered in QCD in the studies of the Regge asymptotics of the scattering 
ampUtudes Ej't3. The eigenvalue problem for the operator Q appears to be much more simple, 
though equivalent, those for the -f/^3/2 Hamiltonian, and a number of results had been obtained in 
the Refs.B'tBllZfllSO for their spectra. In particular, using the different metods it became possible 
to get very accurate decription of the spectrum of the operator Q, and the Hamiltonian H^i2- 
Moreover the obtained results may serve a starting point to analysis of the more interesting from 
the physical point of view case of the A = 1/2 nucleon distribution amplitudes. 

The scale dependence of the A = 1/2 nucleon distribution amplitudes is driven by the 
Hamiltonian 

(5) 



Hi/2 = ^3/2 + y, y 



L7,2 + L-2 



"12 



^23 



The Hamiltonian Hij2 is not integrable and the corresponding eigenproblem cannot be solved 
exactly. Hovewer, the operator V can be treated as a small perturbation to the Hamiltonian 
H^j2 for the most part of the spectrum, except a few low lying levels. The numerical calculation 
gives the spectrum shown on Fig. |l|. The spectra of -^1/2 and H^i2 are very similar in the upper 
part and at the same time two lowest levels of the -^1/2 Hamiltonian appears to be special and 
'dive' considerably below the line of the lowest eigenvalue of -ff3/2- The careful analysis of the 
low energy part of the spectrum | allows to to explain this phenomenon as due to binding of the 
quarks with opposite helicity and forming the scalar diquarkt3. 

The above analysis can be extended to the case of the twist - 3 quark gluon operators, which 
had been paid much attention recently uoeEj. The Q^ evolution of the twist-3 quark gluon 
distribution amplitudes are driven by the Hamiltonians (in the large Nc limit) 

Hs = ^{Jl2 + \)+^{Jl2-\)+^{J2Z + \)+^{J2Z + \)-5, (6) 

Ht = ^{Jl2 + \)+^{Jl2-\)+^{J2Z + \)+^{J2Z-\)-5, (7) 

for the chiral even and chiral odd distributions, respectively, and 5 = 4-0(2) + 3/2. The both 



* It turns out that the resulting problem for the effective Hamiltonian describing dynamics of the low lying 
levels turpi out to be a generalization of the famous Kroning-Penney model of a single particle in a periodic delta 
potential Ej. 



Hamiltonians possess the integrals of motionH which are 

Qs = \Li2,L23j - -L23- -L-i^2^ (8) 

Qt = \Ll2, LI^^ - -LI^ - -L\2- (9) 

The underlying models can be identified as the open nonhoniogeneous noncompact spin chains. 
The theory of the open spin chains is less developed as conmared with those for the closed spin 
chains, and many problems has to be solved yet (see Ref. c3). However many results (such as 
calculation of the asymptotic expansion for conserved charges and energies) are available by the 
elementary methodsoE^I, the full account can be found in Ref.O. 
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